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Abstract 

We have calculated the twodoop strong interaction corrections to the chargino 
pole masses in the DR scheme in the Minimal Supersymmetric Standard Model 
(MSSM) with complex parameters. We have performed a detailed numerical anal- 
ysis for a particular point in the parameter space and found corrections of a few 
tenths of a percent. We provide a computer program which calculates chargino and 
neutralino masses with complex parameters up to 0{aas)- 



1 Introduction 



In the Minimal Supersymmetric Standard Model (MSSM), there are two charginos xt 
and xti which are the fermion mass eigenstates of the supersymmetric partners of the 
W ± and the charged Higgs bosons Hf 2 - Likewise, there are four neutralinos Xi> • • ■ > X°> 
which are the fermion mass eigenstates of the supersymmetric partners of the photon, 
the Z°-boson, and the neutral Higgs bosons H® 2 . Their mass matrices stem from soft 
gaugino breaking terms, spontaneous symmetry breaking in the Higgs sector and in case 
of \i from the super-potential. 

The next generation of future high-energy physics experiments at Tevatron, LHC and a 
future e + e~ linear collider (LLC) will hopefully discover these particles if supersymmetry 
(SUSY) is realized at low energies. Much work has been devoted to the study of the 
physics interplay of experiments at LHC and ILC [Tj. Particularly at a linear collider, 
it will be possible to perform measurements with high precision [U [2J E]. In fact, the 
accuracies of the masses of the lighter SUSY fermions are in the permille region which 
makes the inclusion of higher order corrections indispensible. 

In the framework of the real MSSM important results on quark self-energies were ob- 
tained in [3]-[B]. In [7J[S] the gluino pole mass was calculated to two-loop order. Moreover, 
the MSSM Higgs-sector has been studied in detail, even in the full complex model [5]-[T3]. 

In a previous work [H] we studied loop corrections to the neutralino pole masses 
and found that the relation between the DR-input and the physical observables has to 
be established at least at the two-loop level in order to match experimental precision. 
Following these lines, we calculate in this paper the two-loop 0(aag) corrections to the 
charginos within the MSSM. We conclude that these two-loop corrections are in the mag- 
nitude of the experimental uncertainty and therefore they are relevant e.g. for global fits 
of DR-parameters. 

A new feature of this work is the inclusion of complex parameters in the MSSM. We 
therefore not only study the charginos but also re-analyze the neutralino-masses with 
complex parameters and study in particular the dependence on the phase of the soft 
trilinear breaking parameter A t . 

Generic analytic formulae for SUSY QCD corrections in 0(aas) to fermion pole 
masses in the MSSM were already derived in [TJ. Our calculation is, however, completely 
independent. More precisely, we use semi-automatic Mathematica tools [T5l [16], [TTj for 
the diagram generation and analytic simplifications. 

In the Appendix we briefly describe our C-program Polxino [18] developed for the 
calculation of chargino and neutralino masses with complex parameters up to 0(aas) 
with a convenient SLHA-interface [19| for numerical studies. 

2 Diagrammatics 

In Fig. [1] we show all one-loop diagrams. Similar to the neutralino case we checked our 
analytic one-loop calculation against previous work [20[ [21] in the on-shell scheme and 
found agreement. 

Note that in contrast to the neutralino calculation there are now quark isospin-partners 
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Particle 


Mass 


"LHC" 


"ILC" 


"LHC+ILC" 


X? 

X°2 

xt 

QR 
Ql 
9 


97.7 

183.9 

183.7 

547.2 

564.7 

607.1 


4.8 
4.7 

7-12 
8.7 
8.0 


0.05 
1.2 
0.55 


0.05 
0.08 
0.55 
5-11 
4.9 
6.5 



Table 1: accuracy from experiment at LHC, ILC and LHC © ILC [TJ [21], masses in GeV 

denoted by (q, Q) in the loop which give rise to different tensor reduction formulae. In 
order to get a pure aas correction from (Fig. H|) it is necessary to shorten the 4-squark 
coupling to its QCD part. Diagrams with one- loop counter-term insertions (Fig. [5]) in- 
volve 0(as) mass counter-terms for quarks and squarks as well as coupling constant 
counter-terms stemming from the Yukawa part of the chargino-quark-squark couplings 
and counter-terms to the squark mixing matrix, see e.g. [I]. 

For the evaluation of the amplitudes we adopt the strategy of [2] , that is, we use semi- 
automatic tools [TBI [TBI [TTj in an Mathematica environment and therein auto-create 
Fortran code. 

The renormalization prescription we adopt is the familiar DR -scheme which regulates 
UV-divergencies dimensionally but introduces an unphysical scalar field for any gauge field 
in the theory in order to restore the counting of degrees of freedom in supersymmetry. 
These unphysical mass parameters can be absorbed in the sfermion mass parameters [22J 
(the resulting scheme is called DR') and hence provide a consistency check for the calcu- 
lation of the diagrams containing gluon lines. 

In order to handle infrared divergencies in the individual diagrams we introduce an 
infrared regulating mass parameter and check that the resulting contributions as well as 
the unphysical scalar mass due to the gluon field cancel out in the final result. As the 
main focus of this work is on the numerical analysis we do not reproduce the resulting 
lengthy expressions here. 

Owing to the fact that we split the contributions into self-energy and counter-term 
diagrams, we can quite easily check some of the generic formulae obtained in [7] where 
we find agreement. 

The numerical analysis was performed by implementing Tsil [23] in this Fortran 
program. As in the case of neutralinos we used the usual 't Hooft Feynman i?g=i gauge 
for the gluon field, except for the check of gauge independence. 

3 Numerics 

Our reference scenario used for the numerical analysis is the benchmark point SPSla' 
[21]. The SUSY parameters at Q = 1 TeV are tan/3 = 10, M l = 103.209 GeV, 
M 2 = 193.295 GeV, M 3 = 572.328 GeV, fi = 401.62 GeV, A t = -532.38 GeV, A b = 
-938.91 GeV, M Qa = 470.91 GeV, M &3 = 385.32 GeV and M &3 = 501.37 GeV , for fur- 
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ther details see [23]. The tree-level chargino masses at this point are = 180.9 GeV 
and M-+ = 422.2 GeV. 

Fig. fe] shows the chargino pole masses at SPSla' as functions of tan/?. At the SPSla' 
value tan/3 = 10 we find an absolute two-loop correction 8m~+ = 0.2GeV which is in 
the order of magnitude of the expected experimental uncertainty for this particle, see 
Table [H Therefore, the inclusion of these corrections is mandatory when extracting DR 
-parameters from experiment. 

For Fig. [7J we set the third generation trilinear breaking parameters equal, A 3 = A t = 
Ah. This parameter effects the mixing in the squark sector and therefore enters all the 
two-loop diagrams through the couplings and the sfermion masses. 

In Fig. [8] we assume gauge unification, M± : M 2 : M 3 ~ 1 : 2 : 6. The plot is over M 2 , 
all other values are taken from SPSla'. 

In Fig. [9] we show the one- and two-loop chargino mass shifts as a function of the 
third generation soft SUSY breaking masses Mq = Ma = Mfj 3 = Mp . Again, all other 
parameters are taken from SPSla'. 

Finally, in Fig. [10] we investigate the dependence on <^ t , the complex phase of the 
soft trilinear breaking parameter A t . Extending previous work [T3] we here include the 
neutralinos by taking A t complex. In the first column there are the one-loop corrections 
to the chargino and the neutralino pole masses and in the second the respective two-loop 
corrections. It can be seen that the influence of the phase is quite substantial at the 
two-loop level. 

Fig. HU shows the decrease of the scale dependence when the loop-level of the correc- 
tions is increased. The plots with one- (red) and two-loop (black) masses are zooms of 
the plots to their left where the running tree-level mass is included. The scale dependence 
is reduced considerably when going from the one- to the two-loop level. The remaining 
scaling comes from the uncancelled 0(a 2 ) RGEs. 

4 Conclusions 

We have calculated the chargino pole masses in the MSSM to order O(acts) and performed 
a detailed numerical study. The typical size of the two-loop corrections is comparable to 
the expected eperimental accuracy at future linear colliders and therefore needs to be 
taken into account when analyzing precision experiments. Our analytic expressions agree 
with previous generic results [7J and have been checked thoroughly. Extending previous 
work [13] we also include complex parameters in the neutralino case. Finally, we wrote 
a program [18] for this calculation with an interface to the commonly used SUSY Les 
Houches accord [T9j . 

A Polxino: Short description of the program 

Polxino [18] is a C-program for the evaluation of the pole masses of neutralinos and 
charginos in the MSSM. Up to now it contains the full one-loop level and all SQCD 
effects at the two-loop level. It makes considerable use of the TSIL program [23] and 



4 



uses some of the generic functions of [25]. For convenience we provide a SLHA-interface 
[T9] . We have carefully checked Polxino against the calculation described in this paper 
and in |14j . 



MSSM_DATA member 


description 


gs, gl, g2 


gauge coupling constants 


mue 


/i-parameter 


TB 


tan (3 


VEV 


the square root of the sum 
of squared Higgs VEVs 


MA02 


squared mass of A 


MQU[3], MQD[3], MLE[3] 


quark and lepton-masses 


Ml, M2, M3 


gaugino masses 


MSL[3], MSE[3], MSQ[3], 
MSU[3], MSD[3] 


bilinear soft susy breaking parameters 


Af[4][3] 


trilinear soft susy breaking parameters 



Table 2: basic parameters of the struct MSSM_DATA 



The source code contains the main file main.c which itself contains the code for a 
generic case, the evaluation of aforementioned pole masses at the SPSla' benchmark point. 
Since the program currently uses long double (POL_REAL) and long double com- 
plex (POL_COMPLEX) it is recommended that TSIL is compiled with the compiler 
flag 

-DTSIL_SIZE_LONG 

The central data struct in Polxino is 
MSSM_DATA 

which contains the whole set of MSSM parameters. The subroutine loadSPC fills this 
struct from an SLHA-file and returns 1 if successful. 

loadSPC (" SPheno.spc" , (MSSM_DATA*) ) ; 
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The file SPheno.spc used by the main program was created by SPheno [26] and comes 
with Polxino. If the content of the struct MSSM_DATA is to be changed after the file 
was loaded it must be kept in mind that some of the information stored in the struct is 
derived from a smaller set of basic elements. The basic parameters which may be subject 
to change are given in Table [2] with their respective index ranges in C/C++-notation. 
The first index of Af labels the type of the trilinear breaking parameter according to 
FeynArts. Af[0] denotes the neutrino case and should therefore be kept zero. Af[l], 
Af[2] and Af[3] are the leptonic, up- and down-type set of parameters. All other indices 
of the form [3] label the generation. 
Calling 

mssm_digest((MSSM_DATA*)); 

recalculates all the derived parameters in the struct from the basic ones and re-diagonalizes 
the mixing systems. Thus, it is fairly simple to loop over a parameter that enters in many 
different places, e.g. the //-parameter. The math used for this re-diagonalization was 
extracted from some Fortran subroutines of LoopTools [27] and rewritten in the CD- 
language. 
The command 

mssm_print ( (MSSM_DATA) ) ; 

prints out all the information stored in the struct. Finally, 

getneupole((int) i, (MSSM_DATA), 

(POL_REAL*) Mpolel, (POL_REAL*) Gammal, 
(POL_REAL*) Mpole2, (POL_REAL*) Gamma2); 

getchpole((int) i, (MSSM_DATA) , 

(POL_REAL*) Mpolel, (POL_REAL*) Gammal, 
(POL_REAL*) Mpole2, (POL_REAL*) Gamma2); 

calculate the one- and two-loop pole masses of chargino (neutralino) i and store the one- 
respectivly two-loop results in the variables Mpolel resp.Mpole2. The tree-level and 
one- loop widths are stored in Gammal and Gamma2, respectively. 
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Figure 2: Chargino two-loop selfenergy diagrams with inner gluon line 
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Figure 3: Chargino two-loop selfenergy diagrams with inner gluino line 
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Figure 4: Chargino two-loop selfenergy diagrams with three inner squark lines 
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Figure 5: Chargino two-loop self-energy diagrams with counter-term insertions 
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Figure 6: Absolute chargino mass shifts as functions of tan/3. All other parameters are 
from SPSla'. Left: One-loop mass shifts. Right: Two-Loop mass shifts. In all plots we 
use black for xf an d blue for xt ■ 




Figure 7: Relative chargino mass shifts 8m^+ and 5m-+ as functions of the trilinear 
breaking parameters A 3 = A t = A^. Left: One-loop mass shifts. Right: Two-Loop mass 
shifts. 
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Figure 8: Absolute chargino mass shifts Sm^+ and 5m -+ as functions of the soft gaugino 
breaking mass Mi. In this plot we assume gauge unification, all other parameters are 
from SPSla'. Left: One-loop mass shifts. Right: Two-Loop mass shifts. 
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Figure 9: Absolute chargino mass shifts as functions of Mq, see text. All other parameters 
are from SPSla'. Left: One-loop mass shifts. Right: Two-Loop mass shifts. 
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Figure 10: Absolute neutralino and chargino mass shifts as functions of (f>A t , the complex 
phase of the soft trilinear breaking parameter A t . First column: One- loop corrections. 
Second column: Two-loop corrections. The black line is Xii blue line is \2i red line is x° 
and the green line is %° 
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Figure 11: Scale dependence of chargino pole masses, see text. The blue line is the 
tree-level running DR-mass, red is the one-loop corrected mass and black the two-loop 
corrected mass 
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